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Abstract 



A classical theorem of Stone and von Neumann says that the Schrodinger 
representation is, up to unitary equivalences, the only irreducible represen- 
tation of the Heisenberg group on the Hilbert space of square-integrable 
functions on configuration space. Using the Wigner-Moyal transform we 
construct an irreducible representation of the Heisenberg group on a cer- 
tain Hilbert space of square-integrable functions defined on phase space. 
This allows us to extend the usual Weyl calculus into a phase-space calculus 
and leads us to a quantum mechanics in phase space, equivalent to standard 
quantum mechanics. We also briefly discuss the extension of metaplectic op- 
erators to phase space and the probabilistic interpretation of the solutions 
of the phase space Schrodinger equation 



Contents 



1 Introduction and Motivations! 1 



1.1 General Discussion 2 



1.2 Other possible Schrodinger equations! 3 

1.3 Notations! 6 



2 The Wiener Wave-Packet Transform! 



2.1 Definition and relation with the Wiener- Moval transform! 



2.2 The fundamental property 

2.3 The range of TjA 



3 Phase-Space Wevl Calculus] 13 



3.1 The Schrodinger representation! 13 

3.2 Heisenberg-Wevl operators! 14 



3.3 Extension to phase space! 15 



4 Metaplectic Covariance 18 



4.1 Metaplectic operators! 19 



4.2 Metaplectic group and Wevl calculus! 20 



Schrodinger Equation in Phase Space] 24 

.1 The relationship between ib and J 24 



5.2 Spectral properties! 25 



5.3 The case of quadratic Hamiltoniansl 26 



Interpretation of the phase-space wavefunction m 28 

6.1 Marginal probabilities! 28 

fi.2 The limit h -> fll . ~^ 29 



7 Discussion and Remarks! 



30 



1 Introduction and Motivations 

In a recent Letter ^2] we have shortly discussed and justified the Schrodinger 
equation in phase space 

ih^(x,p, t) = H(x + ihJfe, -ihJt)t>(x,p, t) 

proposed by Torres- Vega and Frederick in j^Hl , and obtained by these 
authors using a generalized version of the Husimi transform. In this paper 
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we set sail to sketch a complete theory for the related equation 
ih^(x,p, t) = H [\x + ih£, \p - ih^j ^{x,p, t). 

where the variables x and p are placed on equal footing. We will see that 
in addition to the greater aesthetic 1 appeal of the latter equation it has the 
advantage of yielding a more straightforward physical interpretations of its 
solutions. 

This paper is reasonably self-contained; we have given rather detailed 
proofs since there are many technicalities which are not always immediately 
obvious. 

1.1 General Discussion 

One of the pillars of non-relativistic quantum mechanics is Schrodinger's 
equation 

dib h 2 

where the operator on the right-hand side is obtained from the Hamiltonian 
function 

H = ^-f + V(r,t) 

by replacing the momentum vector p by the operator — iKV? and letting 
the position vector r stand as it is. But how did Schrodinger arrive at this 
equation? He arrived at it using what the novelist Arthur Koestler called a 
"sleepwalker" argument, elaborating on Hamilton's optical-mechanical anal- 
ogy, and taking several mathematically illegitimate steps (see Jammer ^5] 
or Moore ^jj for a thorough discussion of Schrodinger's argument). In fact 
Schrodinger's equation can be rigorously justified for quadratic or linear po- 
tentials if one uses the theory of the metaplectic group (see our discussion 
in |5] , Chapters 6 and 7), but it cannot be mathematically justified for 
arbitrary Hamiltonian functions; it can only be made plausible by using for- 
mal analogies: this is what is done in all texts on quantum mechanics, and 
Dirac's treatise [I], p. 108-111) is of course not an exception. The gist of 
Schrodinger's argument, recast in modern terms, is the following: a "matter 
wave" consists -as all waves do- of an amplitude and a phase. Consider now 
a particle with initial position vector vq = (x(0), y(0), z(0)). That particle 
moves under the influence of some potential and its position vector becomes 

1 Admittedly, this is a subjective criterion! 
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f(t) = (x(t),y(t), z(t)) at time t. The change of phase of the matter wave 
associated with the particle is then postulated to be the integral 

A$ = - / p- df- Hdt (2) 

calculated along the arc of trajectory T joining the initial point f*o to the 
final point f(tt in space-time; p = {PxiPyiPz) is the momentum vector and 
H = H(f,p,t) the Hamiltonian function. The choice ((2j) for A$ is dictated 
by the fact that it represents the variation in action when the particle moves 
from its initial position to its final position. Now, in most cases of interest 
the initial and final position vectors uniquely determine the initial and final 
momentum vectors if t is sufficiently small, so that HA& can be identified 
with Hamilton's principal function W(fo,f, t) (see 0101); and the latter is 
a solution of Hamilton-Jacobi's equation 

dW 

-Q f + H(r,V f W,t) = 0. (3) 

Schrodinger knew that the properties of the "action form" 

A = p- df - Hdt (4) 

led to this equation, and this was all he needed to describe the time-evolution 
of the phase. We now make an essential remark: the property that A<3? can 
be identified with Hamilton's principal function is intimately related to the 
fact that the action form A is a relative integral invariant. This means that 
if 7 and 7' are two closed curves in the (r , p, t) space encircling the same 
tube of Hamiltonian trajectories, then we have 

p ■ df — Hdt = 6 p-df— Hdt 

1 o // 

(this formula is a consequence of Stoke's theorem and generalizes to an 
arbitrary number of dimensions; see for instance P3E])- 



1.2 Other possible Schrodinger equations 

We now make the following crucial observation, upon which much of this 
paper relies: the action form A is not the only relative integral invariant 
associated to the Hamiltonian H. In fact, for any real scalar A the differential 
form 

A\ = Xp ■ df + (A - l)f • dp - Hdt 
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also satisfies the equality 

& A\ = J> A x 

and is hence also a relative integral invariant. This is immediately checked 
by noting that since 7 is a closed curve we have 



i> p • dr + r ■ dp = i> d(p • r ) = 

J "f J y 



and hence 



(p Xp ■ dr + (A — l)r • dp = (p Xp ■ dr + (1 — X)Xp ■ dr = (p p ■ dr. 

J y 

A particularly neat choice is A = 1/2; it leads to the "symmetrized action" 

A%/2 = 'nip ' dr — r ■ dp) — Hdt (5) 

where the position and momentum variables now play identical roles, up to 
the sign. 

Let us investigate the quantum-mechanical consequences of the choice 
A = 1/2. We consider the very simple situation where the Hamiltonian 
function is linear in the position and momentum variables; more specifically 
we assume that 

Ho =p-r -p - f. 
The solutions of the associated equations of motion 

d -., s _ , d _. . 

— r(t) = r and — = p 

are the functions 

f(t) = r(0) + rot and p(t) = p(0) + pot 

hence the motion is just translation in phase space in the direction of the 
vector (ro,p*o)- An immediate calculation shows that the standard change 
in phase © , expressed in terms of the final position r = f(t) , is 

1 t 2 

A$ = $(r]t) = -(tpo-r-—p -ro)] (6) 
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this function of course trivially satisfies the Hamilton-Jacobi equation for 
Hq . Assuming that the initial wavefunction is ip Q = if) (r) , a straightforward 
calculation shows that the function 



i/j(f, t) = exp 



$(r;t) 



ip (f- tf ) 



(7) 



is a solution of the standard Schrodinger equation 



ih- 



dt 



(-ihf ■V?-p -r)ip = H (f, -ihV r )tp- 



Suppose now that instead of using definition (J2J) for the change in phase we 



use instead the modified action associated with Ai/2- Then 



A$ 1/2 = \ j h(P ' df, ~ r-dp)- Hdt; 



(8) 



integrating and replacing r(0) with r — r§t and p(0) with p — p^t this leads 
to the expression 

®i/2(?,p;t) = -{p- r -po- r) 

which, in addition to time, depends on both r and p; it is thus defined on 
phase space, and not on configuration space as was the case for (jBJ). The 
function &i/2(r,Pl t) does not verify the ordinary Hamilton-Jacobi equation 
(j3J); it however verifies its symmetrized variant 



1/2 



dt 



+ H Q [\f+ V^i /2 , \p- V^ 1/2 ) = 



(9) 



as is checked by a straightforward calculation. This property opens the gates 
to quantum mechanics in phase space: assume again that we have an initial 
wavefunction ip = ipo(f) and set 



ty(r,p, t) = exp 



ip (r- tf ). 



Using Q one finds that 
<9* 



ih- 



dt 



H (if + ihVjj, \p-ihV?) 



(10) 



fill 



there is in fact no reason to assume that the initial wavefunction depends 
only on r; choosing \&o = ^fo{f,p) the same argument shows that the 
function 



^(f,p, t) = exp 



h 



^o(r- tf ,p- tp ) 



(12) 



5 



is a solution of with initial condition *S>q. Observe that the operator 
Hq in the "phase-space Schrodinger equation" is obtained from the 

Hamiltonian function Hq using the phase space quantization rule 

~ d - d 

x — ► X = \x + itiT— , p x — > Px = \px - ih— 

Opx OX 



and similar rules for the y,z variables. The operators X,P X , etc. obey the 
usual canonical commutation relations: 

[x,P x ] = -ih, [Y,Py] = -ih, [z,p~ z } = -ih 

and this suggests that these quantization rules could be consistent with the 
existence of an irreducible representation of the Heisenberg group in phase 
space. This will be proven in Section |3J where we will explicitly construct 
this representation. 

The equation (fTTj) corresponds, as we have seen to the choice A = 1/2 
for the integral invariant A\; any other choice is per se equally good. For 
instance A = 1 corresponds to the standard Schrodinger equation; if we took 
A = we would obtain the phase-space Schrodinger equation 

a* 



ih— = H (f+ ihVf, -ihV r )ty. (13) 



considered by Torres- Vega and Frederick [^Hl HH1 and discussed in ^1 • 

The aesthetic appeal of the Schrodinger equation in phase space in the 
form (|1H) is indisputable, because it reinstates in quantum mechanics the 
symmetry of classical mechanics in its Hamiltonian formulation 

Tt= v '' H ■ ft =" v '* (14) 

in both (P) and (|T1)) the variables x and p are placed, up to a change of sign, 
on the same footing. 



1.3 Notations 

We will work with systems having iV degrees of freedom; we denote the 
position vector of such a system by x = (x%, xn) and its momentum 
vector by p = (pi, ...,pjy). We will also use the collective notation z = (x,p) 
for the generic phase space variable. Configuration space is denoted by 
and phase space by R^. The generalized gradients in x and p are written 

8_ _ ( d d \ d_ _ (d_ d \ 

dx ydxi' ' 8xn J ' dp \dpi' ' dpN J 
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For reasons of notational economy we will write Mu 2 instead of Mu ■ u when 
M is a matrix and u a vector. 

We denote by z A z' the symplectic product of z = (x,p), z' = (x',p'): 

z A z' = p • x' — p' • x 

where the dot • is the usual (Euclidean) scalar product. In matrix notation: 

z A z' = {z') T Jz , J : 



/ 
-I 



where J is the standard symplectic matrix (0 and / are the N x N zero 
and identity matrices). We denote by Sp(N) the symplectic group of the 
(x,p) phase space: it consists of all real 2N x 2N matrices S such that 
Sz A Sz' = z A z'; equivalently 

S T JS = SJS T = J. 

We denote by (•, •) the L 2 -norm of functions on configuration and 
by ((-, •)) that of functions on phase space M. 2N . The corresponding norms 
are denoted by 1 1 • 1 1 and 1 1 1 • 1 1 1 . 

5(M m ) is the Schwartz space of rapidly decreasing functions on R m and 
we denote by F the unitary Fourier transform defined on L 2 (R^ ) by 

mP) = i^nf' 2 I e~^(x)d N x. (15) 

2 The Wigner Wave-Packet Transform 

In what follows <p will be be a rapidly decreasing function normalized to 
unity: 

0€5(M^) , ||0||2 =1. (16) 



2.1 Definition and relation with the Wigner- Moyal trans- 
form 

We associate to (j> the integral operator U$ : L 2 (R^) — > L 2 (R 2Ar ) defined 
by 

U M Z ) = (^h) N/2 e ^ P ' X j e-> x >(x')0(x - x')d N x' (17) 
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and we call Us the "Wigner wave-packet transform" associated with (f>. This 
terminology is justified by the fact that the operator Us is easily expressed 
in terms of the Wigner-Moyal transform 



W(^^)(x,p) = ( 1 L n ) N I e-^y^( x + ly)^ x -ly)d N y (18) 



of the pair (tp,4>) (see 0^1)- I n fact, performing the change of variable 
x' = 7j(x + y) in (JT7|) we get 

/ i \N/2 , 
U^{z) = i^—j 2- N J e-^y^ix + y))HH x ~ vWv 

that is 

u^(z) = {f) N/2 w^Mb). (19) 

Remark 1 A standard -but by no means mandatory- choice is to take for 
<j) the real Gaussian 

*« w = (^)" /4exp ("^ N ' 2 ) ; <20) 

the corresponding operator Us is then (up to an exponential factor) the "co- 
herent state representation" familiar to quantum physicists. 

2.2 The fundamental property 

The interest of the Wigner wave-packet transform Us comes from the fact 
that it is an isometry of L 2 (R^) onto a closed subspace Hs of L 2 (R 2N ) and 
that it takes the operators x and —ihd/dx into the operators x/2 + ihd/dp 
and p/2 — ihd/dx: 

Theorem 2 The Wigner wave-packet transform Us has the following prop- 
erties: (i) Us is an isometry: the Parseval formula 

{{Us^,Us^')) = ^A') (21) 

holds for all Tp,ip' G <S(Mf). In particular U^Us = I on L 2 (R%). (ii) The 
range TLs of Us is closed in L 2 (M. 2N ) (and is hence a Hilbert space), and the 
operator Ps = UsU^ is the orthogonal projection in L 2 (M. 2N ) onto TLs- (Hi) 
The following intertwining relations 

~+iH£)UsiP = UsW) , (l-ih^Us^ = Us{-ih^). (22) 
hold for^ G S(M*f). 
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Proof, (i) Formula (|21j) is an immediate consequence of the property 



((W^AIW^'^'))) = (^f (23) 

of the Wigner-Moyal transform (see e.g. Folland 6 p. 56; beware of the 
fact that Folland uses normalizations different from ours). In fact, taking 
cf) = <j)' we have 

((c/^W)) = (f ) N j w^Mh)wWJ){h)d 2N z 

= (2nh) N (GW,0),W,0))) 

= (^')0^) 

which is formula (|21j) since (j) is normalized. To prove (ii) we note that 
PI = P^ and 

pl = u^up^u; = u;u^ = 

hence P^ is indeed an orthogonal projection. Let us show that its range 
is H<f>; the closedness of TL^ will follow since the range of a projection in 
a Hilbert space always is closed. Since U^U^ = I on L 2 (M^ r ) we have 
Upj^tp = f° r every if; in L 2 ^^) and hence the range of U£ is L 2 (M.^). 
It follows that the range of is that of U^U^ = P^ and we are done. 
(Hi) The verification of the formulae 1)22(1 is purely computational, using 
differentiations and partial integrations; it is therefore left to the reader. ■ 
The intertwining formulae (|22|) show that the Wigner wave-packet trans- 
form takes the usual quantization rules X ^ X j X 

— * ~^~§x lading to the 
standard Schrodinger equation to the phase-space quantization rules 

1 .„ d 1 ,J 

X ^2 X + l % ' P ^2 p - lh dx- ; 

observe that these rules are independent of the choice of (j) and that these 
rules are thus a common features of all the transforms U^. 

2.3 The range of Efy 

One should be aware of the fact that the Hilbert space Ti.^ is smaller than 
L 2 (R 2 Z N ). This is intuitively clear since functions in L 2 (M. 2N ) depend on 
twice as many variables as those in L 2 (W^) of which TL^ is an isometric 
copy. Let us discuss this in some detail. 
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Theorem 3 (i) The range of the Wigner wave-packet transform Lfy asso- 
ciated to the Gaussian \20]) consists of the functions \P G L 2 (M. 2N ) for which 
the conditions 



d . d 

dxj dpj 



for l<j<N (24) 



hold, (ii) For every 4> the range of the Wigner wave-packet transform is 
isometric to Ti.^ . 

Proof, (i) We have = e~2n p ' x V ( j >h where the operator is defined by 



for 1 < j < N. 



It is shown in 22^ that the range of V < f, h consists of all \& G L 2 (M. 2N ) such 
that 

&r l i]) [ OTP (^ lp|a )* (z) 

That the range of is characterized by (|24|l follows by an immediate 
calculation that is left to the reader, (ii) If and are two Wigner 
wave-packet transforms corresponding to the choices 4> 2 then U^U^ is 
an isometry of Ti^, onto Ti^, . m 

The result above leads us to address the following more precise question: 
given ^ G L 2 (R 2 Z N ), can we find <p and t/j in L 2 (R% ) such that ^ = Efy^? We 
are going to see that the answer is no. Intuitively speaking the reason is the 
following: if ^ is too "concentrated" in phase space, it cannot correspond via 
the inverse transform U^ 1 = to a solution of the standard Schrodinger 
equation, because the uncertainty principle would be violated. Let us make 
this precise when the function ^ is a Gaussian. We first make the following 
obvious remark: in view of condition (12411 every Gaussian 



¥o(*) = Aexp^H 2 ^) , AGC 

is in the range of U ( j >h . It turns out that not only does this particular 
Gaussian belong to the range of every Wigner wave-packet transform U^, 
but so does also the compose $0°^ for every S G Sp(N): 

Theorem 4 Let G be a real positive-definite 2N x 2N matrix: G = G T > 0. 
Let *$>g G L 2 (E? Z N ) be the Gaussian: 

9 G (z) = exp (~^GzA . (25) 
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(i) There exist ip,<p & <S(R^) such that U^ip = if and only if G G Sp(N), 
in which case we have 

where (f)^ is the Gaussian A20\) and a an arbitrary complex constant with 
modulus one. (ii) Equivalently, |\Pg| 2 must be the Wigner transform Wip of 
a Gaussian state 



%l>{x) = c (det X) 1/A (ttH) 3N ^ exp 



2ft 



(X + iY)x 1 



(26) 



with \c\ = 1, X and Y real and symmetric, and X > 0. 



Proof. In view of the relation (|19|) between Us and the Wigner-Moyal 
transform ((TU]) is equivalent to 



W(i/>,<t>)(z) = ( ^ 



N/2 



exp 



2 n 

n 



In view of Williamson's symplectic diagonalization theorem [HO] there exists 
S G Sp(N) such that G = S T DS where D is the diagonal matrix 



D 



A 0' 

A 



, A = diag[Ai, X N ] 



the numbers ±£Ai, Xj > 0, being the eigenvalues of JG ; it follows 

that 

/ 2 \ Ar / 2 / 2 
W(^,0)(5'- 1 z) " 



7rft 



exp 



ft 



In view of the metaplectic covariance property of the Wigner-Moyal trans- 
form (see (|47j) in Section 0} there exists a unitary operator 5 : <S(K^) — ► 
5(Mf ) such that 



hence it is no restriction to assume S = I and 



W(^,(j>){z) = ( — 



N/2 



exp 



2 r 



By definition (|18() of the Wigner-Moyal transform this is the same thing as 

2 



(sk)^ / e~>^(x + lyMa - \y)d N y = 2 N exp 



h 



-Dz l 
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that is, in view of the Fourier inversion formula, 



Setting Q = 2A in the generalized Fresnel formula 

N ' 2 j e -iP-y e -^QP 2 d N p = (det Q)-V2 e -^Q- V 
valid for all Q > we thus have 



+ - \y) = 2 N ' 2 (det A)' 1 ' 2 exp 

Setting u = x + \y and v = x — \y this is 



1 



1 



- ( Ax 2 + -Ay 



ip(u)<f>(v) = 2^ (det A) -^x 



exp 



4fr 



((A + A _1 )(m 2 + v 2 ) + 2(A - A^)u • v) 



and this is only possible if there are no terms u ■ v. This condition requires 
that A = A -1 ; since A is positive-definite we must have A = I and hence 
A = I. It follows that 



ip(u)<f>(v) = 2 N ' 2 



exp 



1 

~2h 



(u 2 + v 2 ) 



so that 



V(x)<M0) = fP(0)<f>(x) = 2 N ' 2 exp (-^x 2 



It follows that both ip and 4> are Gaussians of the type 

ij(x)=ii(0)exp(-^-x 2 ) , <f>(x) = <l>(0)ex P (~a?); 



2h 



2h 



since (ft is normalized this requires that (ft = a<ft h with \a\ = 1 and hence 

N/4 



Q Q0 
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bmce we have ^(0)0(0) = 2 N / 2 we must have 



which concludes the proof of part (i) of the theorem. To prove (ii) recall 
from Littlejohn that the Wigner transform of the Gaussian l|2fij) is given 
by the formula 

Wip(z) = exp (- X jGz 2 

where 

_ \X + YX~ X Y YX- 1 ' 
<*- [ X -ly X -l • 

It is immediate to verify that G £ Sp(N) and that G is symmetric posi- 
tive definite. One proves that, conversely, every such G can be put in 
the form above, and which ends the proof of (ii) since the datum of Wi[) 
determines ip up to a complex factor with modulus one. ■ 

3 Phase-Space Weyl Calculus 

Let H/y be the (2N + l)-dimensional Heisenberg group; it is (see e.g. |oll2"ij^ 
the set of all vectors 

(z,t) = (xi,...,XN]Pl,---,PN]t) 

equipped with the multiplicative law 

(z, t)(z', t') = (z + z',t + t' + \z A z). 

3.1 The Schrodinger representation 

The Schrodinger representation of H^v is, by definition, the mapping Ts c h 
which to every (z ,to) in H/v associates the unitary operator Tg c h(zo>*o) on 
L 2 (R^) defined by 

^Sch^O^o)^^) = ex P 

A classical theorem due to Stone and von Neumann (see for instance (fU El] 
for a proof) says that the Schrodinger representation is irreducible (that is, 
no closed subspace of L 2 (R^) other than {0} and L 2 (M^) are invariant by 
^Sch)) and that every irreducible unitary representation of H^r is unitarily 



i 1 

j[-to + Po ■ x - -pq ■ x Q ) i(f(x-xo). (27) 
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equivalent to Ts c h : if T is another irreducible representation of H^r on some 
Hilbert space TL then there exists a unitary operator U from L 2 (W^) to TL 
which is bijective, and such that the following intertwining formula holds: 

[U of Sc -k](z,t) = [f oU](z,t) for all {z,t) in H N . 

Conversely, if U is a unitary operator for which this formula holds, then T 
must be irreducible. We emphasize -heavily!- that in the statement above it 
is nowhere assumed that TL must be L 2 (IR^ r ); it can a priori be any Hilbert 
space, and in particular it can (and will be!) any of the spaces TL^ defined in 
Theorem |21 We will come back to this point in a while, but let us first recall 
how one passes from the Heisenberg group to the Weyl pseudo-differential 
calculus. In Weyl calculus one associates to a "symbol" A = A(x,p) an 
operator A on 5(1^) by the formula 

M{x) = (^- h f J [ e^^A{\{x + y),p)^(y)d N yd N p. (28) 

This formula makes perfectly sense if for instance A E 5(M 2Ar ), and one eas- 
ily verifies that the 
quantization rules 



ily verifies that the "Weyl correspondence" A A leads to the standard 



x ^X = x , P Sp = -4. (29) 

ox 

For more general symbols the double integral must be interpreted in 
some particular way. For instance, if A belongs to the standard symbol 
class S™ S (M.1 N ) with < 5 < p < 1 that is, if for every compact K C 
and all multi-indices a, (3 G there exists Ck,o,/3 such that 



\dSdgA(x,p)\<C Kta ,p(l + \p\ 



\m-p\a\+6\/3\ 



then Q28|) should be viewed as an "oscillatory integral". There are however 
other possible ways to interpret this formula and make it rigorous; we refer 
to jSlinil^ for detailed discussions. (In particular it is shown in jSJ that if 
A is a Hilbert-Schmidt operator if and only if A G I? (M^ ) ) . 

3.2 Heisenberg— Weyl operators 

There is another very useful way of writing Weyl operators, and this will 
lead us to Weyl calculus in phase space. Setting to = in formula 1)27(1 one 
obtains the so-called Heisenberg-Weyl operators Tg c h(zo) : 



Tsch(zo)tp(x) = exp 



jiPo -x - -p ■ Xo) 



i>(x-x ). (30) 
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It is easy to show, using Fourier analysis, that we can write the operator 
(J2SJ) in the form 

WW = {?k) N j (^A)(z )T Sch (z )i>(x)d 2N z (31) 

provided that J- a A, the "symplectic Fourier transform" of A, exists. The 
latter is defined in analogy with the ordinary Fourier transform on R^ by 

^A(z) = (^nf J e-i^'A(z>)d™z>. (32) 

The conditions of existence of J- a A are actually the same as for the usual 
Fourier transform on L 2 (M. 2N ) to which it reduces replacing z by — Jz. No- 
tice that J- a is an involution: T 2 = /. 

Remark 5 It is often convenient to write formula h'Jl)) more economically 
as 

A=(^n) N f {^A)(z)f Sch (z)d 2N z (33) 

where the right-hand-side is interpreted as a "Bochner integral", i.e. an 
integral with value in a Banach space. 



3.3 Extension to phase space 

We now observe that when a Weyl operator is written in the form 1)31 [I or 
(|33j) it literally begs to be extended to phase space! In fact, we can make 
the Heisenberg-Weyl operators act on functions \P in L 2 (M. 2N ) by replacing 
definition (|3()j) by its obvious extension 

T Sc h(zo)^{z) = exp 

and thereafter define the action of A on ^ £ 5(R;^) by the formula 

M*) = (^nf j(^A)(z )T Sch (z Wz)d 2N z . 

This (perfectly honest) choice would lead, using the method we will 
outline in Sectional to the Torres- Vega equation 

ih— =H[x + ih—, -ih— * 
at \ op ox J 



^(Po -x - -Po ■ x ) 



V(Z - Zq) 
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(see |281 I29| ) which we discussed in 12 . Since we prefer, for reasons of 
symplectic covariance, a more symmetric phase-space Schrodinger equation 
in which x and p are placed on equal footing, we will replace To, c h(zo) with 
the operator T p h(zo) given by 

fph(;Zb)¥(z) = exp (—^z A 2hj <K(z - zq) (34) 

(the subscript u ph" standing for phase space) and then define the phase- 
space Weyl operator associated to A by the formula 

= ilk) N J (^A)(z )f ph (z )y(z)d 2N z . (35) 

This operator Ap^ acts continuously on <S(R^) provided that A is a bona fide 
symbol and can hence be extended to L 2 (M. 2 ) by continuity. In accordance 
with the convention in Remark [5] we will often write for short 

4* =(lk) N J (F«A)(z)T ph (z)d 2N z (36) 

where the right-hand side is again viewed as a Bochner integral. 

Observe that the operators T p h satisfy the same commutation relation 
as the usual Weyl-Heisenberg operators: 

T p h(zi)f ph (z ) = exp {^-jZq A zA f ph (z )fph(2i) (37) 



and we have 



T p h{zo)T p h{zi) = exp (j^zo A Zi \ f ph (z + zi). (38) 

These properties suggest that we define the phase-space representation 
f ph of Hiv in analogy with Cfl by setting for $ G L 2 (R 2 Z N ) 

T ph (z ,t )y{z) = ei to T ph (tz )t>(z). (39) 

Clearly T p h(^o 5 *o) is a unitary operator in L 2 (R 2N ); moreover a straightfor- 
ward calculation shows that 

Tph(z ,to)f ph (zi,ti) = f ph (z + zi,t +h + \zq A zi) (40) 
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hence T p h is indeed a representation of the Heisenberg group in L 2 (M. 2N ). 
We claim that the following diagram is commutative for every Wigner wave- 
packet transform Us- 



Tsch i 
L 2 (R?) 



L 2 (R 2 Z N ) 
L 2 (R 2 Z N ). 



More precisely: 



Theorem 6 Let be an arbitrary Wigner wave-packet transform. (i)We 
have 

T P h{zQ,ta)U<i> = U^Tschizo, t ) (41) 

hence the representation T p h is unitarily equivalent to the Schrddinger repre- 
sentation and thus an irreducible representation o/Htv on each of the Hilbert 
spaces Tifj). (it) The following intertwining formula holds for every operator 

Aph- 

AphUtf, = U^Asch- (42) 

Proof. Proof of (i). It suffices to prove formula 1)41(1 for to = 0, that is 

Tph^o)^ = ^T Sch (zo)- (43) 

Let us write the operator in the form = e^ v ' x W^ where the operator 
W(f, is thus defined by 

= (^n) N/2 

We have, by definition of T p h(zo) 



e-tr*tf,(x - x')^(x')d N x'. 



(44) 



Tph(zo)U^ip(z) = exp 



2h 



(z A Zq + (p - p ) ■ (x - Xq) 



W^(z - zo) 



— {-2p- x +p -xo+p-x) 
In 



W^{z - zo) 



= exp 

and, by definition of W^ip, 
W^(z - z ) = {^n) N/2 J e-fc-nWfix -x 1 - xoMx')d N x' 

= {wn) N/2 ^ (P ~ Po) ' X ° [ e~Ti(P~P°> x "-$(x - x")i>{x")d N x'' 
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where we have set x" = x' + xq- The overall exponential in T p h(zo)U ( f > ip(z) 
is thus 



u\ = exp 



Similarly, 



7rr(-Po • xo + P ■ x - 2p ■ x" + 2p • x") 
2fi 



e -W x "<t>{x - x")eTM*"-2P^)^( x " - XQ )d N x" 



yielding the overall exponential 



U2 = exp 



{ 2 P-x-p- x" + po ■ x" - -po ■ x 



Ul 



which proves (|43|) . The irreducibility statement follows from Stone-von 
Neumann's theorem. Let us prove formula (|42|) . In view of formula (|4.3|) we 
have 

W = (^h) N J (^A)(z )T ph (z )(U^)(z)d 2N z 
= iwn) N J (^A)(z )C/ (/) (fs ch (zo)^)(z)d 2Jv z 
= i^nfUf (/ (F a A)(z )T Sch (z )^(z)d 2N z 

= ^(IscWOO) 

(the passage from the second equality to the third is legitimated by the fact 
that U<f, is both linear and continuous). ■ 



4 Metaplectic Covariance 

Since Sp(N) is the symmetry group for the usual CCR (canonical commu- 
tation relations) 

[Xj,P k ] = ihS jk 

for Xj = Xj, Pk = —ihd/dxk), the uniqueness of these relations implies that 
for each S in Sp(N) there must be some associated unitary operator linking 
them to quantization. These operators are the metaplectic operators. 
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4.1 Metaplectic operators 

Let us recall how the metaplectic operators are defined (for details and proofs 
see for instancejOllHliniE])- Assume that S is a free symplectic matrix, that 
is S G Sp(N) and 



S 



A B 
C D 



with det B ^0. 



To every such S one associates the operators ±Sw,m defined by the formula 



S w ,mip(x) 



i \ n/2 im 

^ih) y/\detB\ 



'} w ^ x '^(x')d N x'; (45) 



here W is Hamilton's characteristic function familiar from classical mechan- 
ics (see for instance 0111): 



W(x, x') = ^DB^x 2 - B~ l x ■ x' + ^B^Ax' 2 , 



(46) 



and m is an integer ("Maslov index") corresponding to a choice of argdet B. 
The operators Sw,m are a sort of generalized Fourier transform, and it is 
not difficult to check that they are unitary. In addition the inverse of S\y,m 
is given by 



S 



W.m 



S w* 



, W*(x,x) = -W(x',x) , m*=N-m 



hence the Sw,m generate a group: this group is the metaplectic group Mp(N) 
(see de Gosson 9 for a complete discussion of the properties of Mp(N) and 
of the associated Maslov indices) . If we choose for S the standard symplectic 
matrix 

" I 
-I 



J 



the quadratic form (|Ht|) reduces to W (x, x') = —x ■ x'\ choosing arg det B 
arg det 1 = the corresponding metaplectic operator is 



Jip(x) 



1 xAT/2 

2nih J 



e-i x - x '^(x')d N x' = i- N l 2 F^(x') 



where F thus is the usual unitary Fourier transform 151) . 

The Wigner-Moyal transform enjoys the following metaplectic covari- 
ance property: for every S S Mp(N) with projection 8 E Sp(N) we have 



(47) 
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Since the Wigner wave-packet transform is defined in terms of W(ip, 0) by 
formula 1)19|) it follows that 

j- \ N/2 _ 



) W^S-^iS-'iz)) 



and hence 



U^SiP) = {Ufr.il>) ° , <fa = S-i<f>. (48) 



4.2 Metaplectic group and Weyl calculus 

In ^3] we have shown, following an idea of Mehlig and Wilkinson JJj, that 
the metaplectic group is generated by operators of the type 

/ i \N/2 MS) r 

S * (I, = (^) ^S-J ^"^ 1 "^^" <49) 

where det(S* — /) / 0, Ms is the symplectic Cayley transform of S: 

M s = ^J(S + I)(S-I)-\ 

and v(S) is the Conley-Zehnder index (modulo 4) of a path joining the 
identity to / in Sp(N) (see for instance Muratore-Ginnaneschi |2J for a 
discussion of this index). For instance, if S = Sw,m then 

^ / i \ N / 2 ,-m-Inert /• 

S^,) = (_) 7mrWl J .-««Il U .(»W.)d»» (50) 

where Inert the number of negative eigenvalues of the Hessian matrix 
of W. Formulae (|49|) and (|5U|) are thus the Weyl representations of the 
metaplectic operators S and Sw,m)- They allow us to define phase-space 
metaplectic operators S p h in the following way: if det(5 — I) ^ we set 

/ 1 \ N / 2 MS) r . 

S « MZ) = (SS) V |de t (5-/)| / CW,^ (M) 



20 



the operators S p h are in one-to-one correspondence with the metaplectic 
operators S and thus generate a group which we denote by Mp p ^(N) (the 
"phase space metaplectic group"). In following lemma we give alternative 
descriptions of the operators (|49|) in terms of the operators T p h: 

Lemma 7 Let S p h £ Mp p h(N) have projection S £ Mp(N) such that 
det(S-I) + 0. We have 

S P k={^) N i u{S) V\det(S-I)\ J e-i SzAz f ph ((S-I)z)d 2N z (52) 

and 

V = {^) N i v(S) V\det{S-I)\ J T ph {Sz)T ph {-z)d 2N z. (53) 
Proof. It is mutatis mutandis the proof of Lemma 1 in JH] ) : we have 

§ J(S + I)(S- i)- 1 = \J + J(S - 1)- 1 

hence, in view of the antisymmetry of J, 

M s z -z = J{S- iy x z -z={S- iy x z A z 

Performing the change of variables z i — > (S — I)z we can rewrite the integral 
in the right-hand side of ((oTl as 

J e ms^)f(z)d 2N z = y/\det(S-I)\ J e^ s - I >f ph ((S-I)z)d 2N z 

= V|det(5-/)| J e- L 2 SzAz f ph ((S-I)z)d 2N z 

hence (|52|). Taking into account formula ([rSSjl for the product of two meta- 
plectic operators T p h(zo) and T p h(^i) we get 

T((5 - I)z) = e^ SzAz f ph (Sz)f ph (-z) 

and formula (|53j) follows. ■ 

This result will allows us to show in a simple way that the well-known 
"metaplectic covariance" relation 

A^Ssch = S-^schS (54) 

valid for any 5 E Mp(N) with projection 5 € Sp(N) extends to the phase- 
space Weyl operators A p ^ provided one replaces Mp(N) with Mp p h(N). 
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Theorem 8 Let S be a symplectic matrix and S p h any of the two opera- 
tors in Mp p h(N) associated with S. The following phase-space metaplectic 
covariance formulae hold: 

Sph,Tph(zo)S~ h = T ph (Sz) , A o Sph = Sp^AphSph- (55) 

Proof. To prove the first formula (|55|) it is sufficient to assume that det(<S — 
I) 7^ and that .Sph is thus given by formula (|5Tj) : since such operators 
generate Mp p h(N). Let us thus prove that 

f ph (Sz )S ph = S ph f ph (z ) if det(5 -1)^0. (56) 
Using either formula (|53[) in Lemma [7| above and setting 

Cs=(^)V(*y|det(S-/)| 

we have 



2N Z 



and 



Setting 



and 



T p h(5^o)^ P h = C s J f ph (Sz )f ph (Sz)f ph (-z)d 
S p bT p h(zo) = C s J f ph (Sz)f ph (-z)f ph (z )d 2N z. 
A(z ) = J T ph (Sz )f ph (Sz)T ph (-z)d 2N z 



B(zo) = J f ph (Sz)f ph (-z)f ph (z )d 2N z 
we have, by repeated use of 

A(z ) = J e^ z ^f ph (Sz + (S-I)z)d 2N z 
B(z ) = Je^ z ^T ph (z + (S-I)z)d 2 

where the phases 3>i and $2 are given by 

$i(z,z Q ) = z A z - S(z + z ) A z 
®2(z,zq) = -Sz A z + (S - I)z A z . 
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Performing the change of variables z' = z + zq in the integral defining A(zq) 
we get 

A{z ) = J e^ l(z '- zo ^T ph {zo + (S- I)z')d 2N z' 

and 

$1(2' - Zq, Zq) = Z A (z' - Zq) - Sz A (z' - Zq) 

= (S — I)z' A Zq — Sz' A z' 

= ®2(z',Zq) 

hence (|56[), The second formula (|55|1 easily follows from the first: noting 
that the symplectic Fourier transform (|32j) satisfies 

F*[AoS\{z) = (^n) N Je-i^'A(Sz')d™z> 

= (^f Je-i Sz ° A *'A(z')d 2N z> 

= FaA(Sz) 

we have 

= {ikf J F a A(Sz)f ph (z)d 2N z 
= (wnf J^A(z)f ph (S- 1 z)d 2N z 
= (wn) N J^A(z)S^f ph (z)S ph d m z 

which concludes the proof. ■ 

It can be shown, adapting the proof of a classical result of Shale [2E] (see 
Wong [HJ, Chapter 30, for a proof) that the metaplectic covariance formula 

A O Sph = 'S'ph^-ph'S'ph 

actually characterizes the phase-space Weyl operators A p ^. That is, any op- 
erator satisfying this relation for all operators S p ^ G Mp p ^(N) is necessarily 
of the type 

A* = {£*ff (^A)(z)f ph (z)d 2N z. 

For example, if H is the Hamiltonian function of the one-dimensional 
harmonic oscillator put in normal form 

H = ^(p 2 + x 2 ) (57) 
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we get 

^ P h = -— V2- ?T zAV, + -|z| 2 (58) 
where V z is the gradient operator in (x,p). 

5 Schrodinger Equation in Phase Space 

We now have all the machinery needed to justify and stud the Schrodinger 
equation in phase space. 

5.1 The relationship between ip and ^ 

The following consequence of theorem links standard "configuration space" 
quantum mechanics to phase-space quantum mechanics via the Wigner wave- 
packet transform and the extended Heisenberg group studied in the previous 
sections. For clarity we denote by ^4sch the usual Weyl operator associated 
by to an observable A. 



Corollary 9 LetU^, <ft S S(M.^), be an arbitrary Wigner wave-packet trans- 
form, (i) If tp = ip(x,t) is a solution of the usual Schrddinger's equation 

= HschV 

at 

then m = (C/^VX 2 ^) is a solution of the phase-space Schrodinger equation 

dm 

ih— = H p fi$>. (59) 

(ii) Assume that m is a solution of this equation and that = ^("jO) 
belongs to the range TCa of Ua. Then m(-,t) 6 Ha for every t for which m 
is defined. 

Proof. Since time-derivatives obviously commute with Ua we have, using 
dm 

ih— = Efy(-HsdW0 = HphtU+ip) = H ph m 

hence (i). Statement (ii) follows. ■ 

The result above leads to the following interesting questions: since the 
solutions of the phase-space Schrodinger equation l|59*|) exist independently 
of the choice of any isometry Ua, what is the difference in the physical 
interpretations of the corresponding configuration-space wavefunctions ip = 
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U^S> and ip' = U^,^? The answer is that there is no difference at all provided 
that <j) and <p are not orthogonal: 

Theorem 10 Let ^ be a solution of the phase space Schrddinger equation 
\5 ( J\) with initial condition and define functions tp 1 and tp 2 ^ n ^(M-x) by 

* = c/^Vi = u^ 2 . 

We assume that $o G Ufa H Tifa ; 

(i) We have G 7^ n /or a// i 

(m,) 7/ (^1,^2) = i/ien Vi anc ^ V'2 are orthogonal quantum states: 

(V!,V 2 ) = o. 

Proof. Property ^ follows from Theorem E3fm). Let us prove (ii). In view 
of formula (|2*5|) we have 

that is H^IH 2 = A(V'i,V ; 2) with A = (<^i,0 2 )- The assertion follows. ■ 
5.2 Spectral properties 

The operators A p ^ defined by Q35|) enjoy the same property which makes the 
main appeal of ordinary Weyl operators, namely that they are self-adjoint 
if and only if their symbols is real. 

Theorem 11 Let A p h and A$ c h be the operators associated to a symbol A. 
We assume that the symplectic Fourier transform T a A is defined, (i) The 
operator A p h is self-adjoint in L 2 (R 2N ) if and only if A = A; (ii) Every 
eigenvalue of As c h is also an eigenvalue of A v ^. 

Proof, (i) By definition of A p ^ and T p h we have 

= i^hf j^A(z )e-^ Az ^(z-z )d 2N z 
= {^f J^A(z-z')e^ Az '*(z')d 2N z' 

hence the kernel of the operator A p ^ is 

K(z,z') = (^ E ) N e^'^A(z-z>). 
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In view of the standard theory of integral operators A p ^ is self-adjoint if and 
only if K{z,z') = K(z',z); in view of the antisymmetry of the symplectic 
product we have 

~W^j = {^n) N e^ zAz '^A(z' - z) 
hence our claim since by definition (|32l) of the symplectic Fourier transform 

T^MP ~ z) = {^n) N J e-^-^'Ajz^d^z" = TjA{z - z'). 

(ii) Assume that A^tp = choosing <p G 5(R^) we have, using the 
intertwining formula Q42JI . 

hence A is an eigenvalue of A p ^. m 

Notice that there is no reason for an arbitrary eigenvalue of ^4 p h to be 
an eigenvalue of A^ch', this is only the case if the corresponding eigenvector 
belongs to the range of a Wigner wave-packet transform. 

5.3 The case of quadratic Hamiltonians 

There is an interesting application of the theory of the metaplectic group 
outlined in Section 0] to Schrodinger's equation in phase space. Assume 
that H is a a quadratic Hamiltonian (for instance the harmonic oscillator 
Hamiltonian); the flow determined by the associated Hamilton equations is 
linear and consists of symplectic matrices St- Letting time vary, thus obtain 
a curve 1 1 — ► St in the symplectic group Sp(N) passing through the identity 
/ at time t = 0; following general principles to that curve we can associate 
(in a unique way) a curve t i — > St of metaplectic operators. Let now 
■00 = V"o( x ) be some square integrable function and set tp(x,t) = Stip (x). 
Then tp is just the solution of the standard Schrodinger's equation 

ih^ = HiP , V(i = 0)=Vo (60) 

associated to the quadratic Hamiltonian function H. (Equivalently, St is 
just the propagator for (JHOJ).) This observation allows us to solve explicitly 
the phase-space Schrodinger equation for any such H. Here is how. Since 
the wave-packet transform U automatically takes the solution ip of H6U|) to 
a solution of the phase-space Schrodinger equation 
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we have 



*{z,t) = (S t ) ph V(z,0). 



Assume now that the symplectic matrix St is free and det(S^ — /) 7^ 0. Then, 
by (ED, 



1 \ N/2 -mit) -Inert W xx (t) r . „ 

' * ' 1 e^ M ^f ph (z )*(z,0)d 2N z 



j\tet{s t -i)\ 



(61) 

where m(t), W xx (t), and Ms(t) correspond to St- If t is such that St is 
not free, or det(<% — I) = 0, then it suffices to write the propagator St as 
the product of two operators (|50[): note however that such values of f are 
exceptional, and that the solution (|61|) can be extended by taking the limit 
near such t provided that takes some care in calculating the Maslov indices. 

Let us illustrate this when H is the harmonic oscillator Hamiltonian 
function (|57j) . The one-parameter group (St) is in this case given by 

cos cot sin cot 
— sin cot cos cot 

and the Hamilton principal function by 



W(x,x';t) 



1 



2 sin cot 

A straightforward calculation yields 



((x + x' ) cos cot — 2xx'). 



M s (t) 



—2 cos ut+2 









-2 cos uit+2 



cot 



cot(f ) 



and 



det(5t - I) = 2(1 - coscot) = 4sin 2 (f ); 



moreover 



W xx (t) = -tan(f )• 
Insertion in formula (|filj) yields the explicit solution 



*(z,t) 



2 2-7r/isin 



•^|l/2 



cxp 



i . 2 2\ f U ^\ 

^:(x +p )cot(y / 



T ph (zo)*(z,0)d 2 zo 



with 



if < t < I 
-2 if - ? < t < 0. 
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6 Interpretation of the phase-space wavefunction 



Let us shortly discuss the probabilistic interpretation of the solutions \& of 
the phase-space Schrodinger equation 

we will in particular elucidate the role played by <j>. 
6.1 Marginal probabilities 

Let ip be in L 2 (B^); if ip is normalized: \\ip\\ = 1 then so is ^ = U^ip in view 
of the Parseval formula (|21|) : IH^IH = 1. It follows that l^l 2 is a probability 
density in phase space. It turns out that by an appropriate choice of <f> the 
marginal probabilities can be chosen arbitrarily close to \ip\ 2 and l-F^I 2 : 

Theorem 12 Let tp G L 2 (R%) and set * = U$ip. (i) We have 

\*{x,p)\ 2 d N p={\<P\ 2 *\iP\ 2 ){x) (62) 

W(x,p)\ 2 d N x = (\F(f>\ 2 * \FiP\ 2 )(p). (63) 



/ 



(ii) Let (A),^ = (Aschrfiip) be the mathematical expectation of the symbol A 
in the normalized quantum state ip. We have 

(A)j = ((A ph %*)) , 9 = U+il>. (64) 

Proof. We have, by definition of 

^ (Z)|2 = (i) / / e_ ^' (X '~ X '' ) ^(^- x 0^(^-^ / )V'(^)?(^ / )d A 'x / d Jv x // . 
Since we have, by the Fourier inversion formula, 

J e -y< x '- x "^d N p = (2irh) N 5(x> - x") 

it follows that 

\^(z)\ 2 d N p = [ff 8{x' - x")\<P(x - x')\ 2 \4>{x')\ 2 d N x'd N x" 

5{x' - x")d N x" \<j){x - x')\ 2 \iP(x')\ 2 d N x' 
\<P{x-x')\ 2 \iP{x')\ 2 d N x' 
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hence formula Q62|) . To prove ()63|) we note that in view of the metaplectic 
covariance formula (|48|1 for the wavepacket transform we have 

U jj,(Ji>)( x ,p) = u,ptp(-p,x) 

where J = i~ N / 2 F is the metaplectic Fourier transform. It follows that 

U F4> {F$){x,p) =r N U^(-p,x) 
and hence changing (— p, x) into (x,p): 

U^{x,p)=i N U F ^{Fi,){p, -x). 
and hence, using (|62|). 

y \*(x,p)\ 2 d N x = J \U F(j> {F^)(p,-x)\ 2 d N x 

\U F4> {FiP)(p,x)\ 2 d N x 

= (\F<p\ 2 *\FiP\ 2 )(p) 

which concludes the proof of (|63|) . To prove (|64|) it suffices to note that in 
view of the intertwining formula (|42|) and the fact that U± = U^ 1 we have 

((Aph*,*)) = ((A ph U^,U^)) 

= (u;A ph u^,4>) 

m 

The result above shows that the marginal probabilities of |^| 2 are just 
the traditional position and momentum probability densities \i/j\ 2 and l-F^I 2 
"smoothed out" by convoluting them with |<^>| 2 and |-F<?!>| 2 respectively. 

6.2 The limit H -> 

Assume now that we choose for <f) the Gaussian H2U|) : 

2 



2h' xl 



The Fourier transform of d> is identical to 
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hence, setting ^ = TJj, ip, and observing that |0^| 2 — > 6 when h — > 0: 

J \^ h {x, P )\ 2 d N p = (h 2 * i^i 2 )(x) ^ m X )\ 2 

J \* H {x,p)\ 2 d N x = (\F^\ 2 * \^ h \ 2 )(p) ^° \F^(p)\ 2 . 

Thus, in the limit h — > the square of the modulus of the phase-space 
wavefunction becomes a true joint probability density for the probability 
densities \ip\ 2 and \Fx/j\ 2 . 

7 Discussion and Remarks 

We have exposed some theoretical background for a mathematical justifica- 
tion of the phase-space Schrodinger equation 

ih^(x,p, t) = H [\x + ih£, \p - ih-^j t). 

The aesthetic appeal of this equation is obvious -at least if one likes the 
Hamiltonian formulation of mechanics. But is this equation useful? While 
the notion of "usefulness" in Science more than often has a relative and 
subjective character, one of the main practical appeal of the phase-space 
Schrodinger equation is that it governs the quantum evolution of both pure 
and mixed states, while the solutions of the usual Schrodinger equation are, 
by definition, only pure states. Another of the advantages of the phase-space 
approach is, as pointed out in JH]) the availability of factorization methods 
for the Hamiltonian, for instance SUSY. From a practical point of view it 
could be held against Schrodinger equations in 2A r -dimensional phase space 
that they are uninteresting because they involve solving a partial differential 
equation in 2iV + 1 variables instead of N + 1 as for the ordinary Schrodinger 
equation. But this is perhaps a somewhat stingy reservation especially in 
times where modern computing techniques allow an efficient processing of 
large strings of independent variables. 

It would perhaps be interesting to make explicit the relationship between 
the theory of Schrodinger equation in phase space we have sketched and 
other approaches to quantum mechanics in phase space, for instance the 
"deformation quantization" of Bayen et al. [H], and whose master equation 
is the "quantum Liouville equation". 

We would like to end this section -and paper!- by discussing a little bit 
the possible physical interpretation of the phase space Schrodinger equation. 
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Recall that we showed in Theorem 0] that a phase-space Gaussian 

*G(z) = exp (-^/) , G = G T >0 

is in the range of any of the Wigner wave-packet transforms Ua if and only 
if G G Sp(N), and that in this case 

W^(z) = \* G (z)\ 2 = exp (-\gz 2 

for some (pure) Gaussian state ip. Let us more generally consider Gaussians 



* M (z) = exp (-^Mz 2 ^j 



where M is a positive-definite symmetric real matrix. One proves that *$>m 
is the Wigner transform W(p) of a (usually mixed) quantum state if and 
only if M _1 + i J is positive-definite and Hermitian: 

{M~ l +iJ)* = M~ l +iJ > 0. (65) 

The probabilistic meaning of this condition is the following: defining as usual 
the covariance matrix of the state p by 

S = -M- 1 
2 

condition (|65|) can be rewritten as 

(S + i^J)* = S + i^J > (66) 

which is equivalent to the uncertainty principle of quantum mechanics (see 
|26l I27j ; we have also discussed this in ^Oj) . For instance, when N = 1 the 
matrix 

Ax 2 A(x,p) 
A(x,p) Ap 2 



E 

satisfies (|66]l if and only if 



Ax 2 Ap 2 > ^k 2 + A(x,p) 

which is the form of the Heisenberg inequality that should be used as soon as 
correlations are present, and not the usual text-book inequality AxAp > ^h. 
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It turns out that conditions (|65|) - Q66JI have a simple topological inter- 
pretation: we have shown in previous work of ours |1(J1 lllj that they are 
equivalent to the third condition: 

The phase-space ball B{y/h) :\z\<h can be embedded into the "Wiener 
ellipsoid" W M ■ Mz 2 < H using symplectic transformations (linear or not). 
Equivalently: the symplectic capacity (or "Gromov width") of Wm is at 
least ttH = ^h, one half of the quantum of action: c(Wm) > \h. 

We have discussed in some detail in [101 lllj how this result allows a 
"coarse graining" of phase space by symplectic quantum cells, which we 
dubbed "quantum blobs" . It appears that it is precisely this coarse-graining 
that prevents Gaussians m with Wigner ellipsoids smaller than a "quantum 
blob" to represent a quantum state. Is this to say that if the Wigner ellip- 
soid of has exactly symplectic capacity \h then f m is a pure state? No, 
because such states are characterized by the fact that the associated Wigner 
ellipsoid is exactly the image of the ball B(sfh) by a symplectic transforma- 
tion since they are described by the inequality Gz 2 = (Sz) 2 < h in view of 
Theorem |11 and there are infinitely many ellipsoids with symplectic capac- 
ity \h which are not the image of B(y/h) by a symplectic transformation. 
However, we have shown in jlUl ITT] that if the ellipsoid Wm '■ Mz 2 < h 
has symplectic capacity \h then one can associate to Wm a unique pure 
Gaussian state. The argument goes as follows: if c(Wm) = \h then if S and 
S' in Sp(N) are such that 



then there exists R E U(N) = Sp(N)nO(2N) such that S = RS' (the proof 
of this property is not entirely trivial) and hence S T S = (S') T S'. It follows 
that the ellipsoid Wq '■ Gz 2 < h, G = S T S, is uniquely determined by Wm 
and that the pure Gaussian state corresponding to 



is canonically associated to the mixed state ^m, which does not in general 
belong to the range of any Wigner wave-packet transform U^. It would 
be interesting to generalize this result to arbitrary functions ^ £ L 2 (M. 2N ) 
by defining, in analogy with the Gaussian "Wigner set" Wqi associ- 

ated with fy. One could then perhaps prove that represents an arbitrary 
(mixed) quantum state provided that W\^ has a symplectic capacity at least 
\h. But enough is enough! We hope to come back to these possibilities in 
the future. 



S(B(Vh)) c W M , S'(B(Vh)) c W M 
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